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Abstract 



Production of S~ and A(1520) in hadronic Z decays has been measured using 
the Delphi detector at Lep. The S~ is directly reconstructed as a charged 
track in the DELPHI microvertex detector and is identified by its S ^ utt de- 
cay leading to a kink between the S~ and vr-track. The reconstruction of the 
A(1520) resonance relies strongly on the particle identification capabilities of 
the barrel Ring Imaging Cherenkov detector and on the ionisation loss measure- 
ment of the TPC Inclusive production spectra are measured for both particles. 
The production rates are measured to be 

{Nj^-fN^""^) = 0.081 ± 0.002 ± 0.010 , 
(iV'A(i520)/iVz^'^) = 0.029 ± 0.005 ± 0.005 . 

The production rate of the A(1520) suggests that a large fraction of the stable 
baryons descend from orbitally excited baryonic states. It is shown that the 
baryon production rates in Z decays follow a universal phenomenological law 
related to isospin, strangeness and mass of the particles. 
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1 Introduction 

The study of baryon production provides an important tool to test models of the 
fragmentation process Beyond the cluster fragmentation model 0, and the string 
model which employs many parameters to describe baryon fragmentation, thermody- 
namical ^ and phenomenological models have appeared recently which success- 

fully describe the overall particle production rates in high energy interactions with very 
few parameters. 

At Lep it has been shown that a large fraction of the observed stable mesons stem 
from decays of scalar and tensor mesons with angular momentum. For baryons this is not 
the baryon resonances, especially those with orbital excitations, typically have 

a large decay width and complicated decay modes. Hence these states are difficult to 
access experimentally in a multihadronic environment. In any case, it is still a question 
of basic importance as to how far baryon production leads to excited baryonic states. 

So far the only orbitally excited baryonic state measured in e^e~ annihilation is 
the A(1520) 0,^]. This paper provides further data on A(1520) and S" production[| 
in hadronic Z decays. For production at the Z so far only two measurements are 
available |lOJTl|, one of them being sensitive only to the sum of S~ and E"*" states |10 . 



The data used throughout this paper were collected by the Delphi detector in 1994 
and 1995. In these data taking periods the Delphi micro Vertex (VD) and Ring Imaging 
Cherenkov (RICH) detectors were optimally set up and functioning for the analyses 
presented. 

This paper is organised as follows. Section gives a brief overview on the detector, 
experimental procedures used to select tracks and hadronic events as well as on the 
specific experimental procedures and corrections used for S~ and A(1520) reconstruction. 
Section (||) presents the results for S~ and A(1520) production and the corresponding 
systematic errors. These results are compared to the expectation of fragmentation models 
and a general phenomenological law of baryon production in Z decays is deduced. Finally 
we conclude in Section (|). 



2 The Experimental Procedure and Event Selection 



The Delphi detector is described in detail in []12|. The present analysis relies on 



information provided by the central tracking detectors and the barrel RICH: 

• The micro Vertex Detector (YD) consists of three layers of silicon strip detectors 
at radii of 6.3, 9.0 and 10.9 cm. Rcf) coordinates^ in the plane perpendicular to the 
beam are measured in all three layers. The first and third layers also provide z 
information. The polar angle [9) coverage for a particle passing all three layers is 
from 44° to 136°. The single point resolution has been estimated from real data 
to be about 8 /im in Rcf) and (for charged particles crossing perpendicular to the 
module) about 9 /im in z. 

• The Inner Detector (ID) consists of an inner drift chamber with jet chamber 
geometry and 5 cylindrical MWPC (in 1995 straw tube) layers. The jet chamber. 



'^The antiparticles are always implicitly included. 

■^In the standard DELPHI coordinate system, the z axis is along the electron direction, the x axis points towards the 
centre of LEP, and the y axis points upwards. The p olar ang le to the z axis is denoted by 6, and the azimuthal angle 
around the z axis by 0; the radial coordinate is R = y x'^ + y'^ . 
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between 12 and 23 cm in R and 23° and 157° (15°-165° for 1995) in 9, consists of 24 
azimuthal sectors, each providing up to 24 R(f) points. 

• The Time Projection Chamber (TPC) is the main tracking device of Delphi. 
It provides up to 16 space points per particle trajectory for radii between 40 and 
110 cm. The precision on the track elements is about 150 /xm in Rcj) and about 600 
Hm in z. A measurement of the energy loss dE/dx of a track is provided with a 
resolution of about 6.5%. 

• The Outer Detector (OD) is a 4.7 m long set of 5 layers of drift tubes situated at 
2 m radius to the beam which provides precise spatial information in Rcf). 

• The Barrel Ring Imaging Cherenkov Counter (BRICH) is the main Delphi 
detector devoted to charged particle identification. It is subdivided into two halves 
[z ^ 0) and provides particle identification using Cherenkov radiation produced 
in a liquid or a gas radiator. This radiation, after appropriate focusing, is trans- 
formed into photoelectrons in a TPC-like drift structure and the Cherenkov angles 
of the track in both media are determined. The BRICH detector provides particle 
identification in the momentum range 0.7 to 45 GeV/c. 

An event was selected as a multihadronic event if the following requirements were satisfied: 

• There were at least five well measured charged particles in the event, each with 
momentum larger than 300 MeV/c (400 MeV/c for the S~ analysis) and a polar 
angle in the range 20° < 9 < 160°. 

• The total reconstructed energy of these charged tracks had to be larger than 12% of 
the centre-of-mass energy. 

• The total energy of the charged particles in each detector hemisphere (defined by 
the plane perpendicular to the beam axis) had to exceed 3% of the centre-of-mass 
energy. 

• The tracking devices and, in the case of the A(1520) analysis also the BRICH, were 
fully operational. 

After these cuts, about 1.3 million events remained for the 1994 period and 0.6 million 
events for the 1995 run around the Z° pole. The A(1520) analysis is based on both years 
data, the S~ analysis only on the 1994 data. 

To study the influence of cuts, inefficiencies and resolution as well as particle re- 
interactions in the detector, a large set of simulated Z — qq events has been used. This 
simulated sample has been generated using the Jetset 7.3 model |^ with the parton 
shower option. The model parameters were taken from earlier QCD studies |]TB[. The 



initial event simulation was followed by a detailed detector simulation |T^. For the 
A(1520) study a specific set of 10000 of these events has also been produced with at least 
one A(1520) — > pK~ per event. 



2.1 Reconstruction 

The charged S hyperons decay through the weak interaction according to 

PttO (^ 52%), 

n7r+ (^48%), (1) 
nvr- (^ 100%). 

The reconstruction of the decay rni^ is based on the large fiight length (cr = 

4.43 cm for S~ and cr = 2.4 cm for S"^). It allows a determination of the track param- 
eters for the S^, if there are at least three hits in the microvertex detector. The decay 
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S — ^ nvr is then reconstructed by finding the kink between the S and the pion which 
is normally well measured in the other tracking chambers, especially the TPC Thus, the 
detection of the neutron is not necessary. The tools needed for the analysis are described 
in more detail below. 

2.1.1 Track reconstruction with the Microvertex Detector 

The geometry of the 1994 Delphi microvertex detector allows the determination of 
the track parameters of a charged particle in two or three dimensions using information 
from this detector alone. These tracks, called 'VD tracks' in the following, can arise from: 

• Low-momentum charged particles {p < 50 MeV/c) not reaching the TPC due to the 
bending in the magnetic field. 

• Tracks at the borders of the TPC modules not efficiently reconstructed by the stan- 
dard tracking algorithm. 

• Charged particles interacting with the detector material outside the VD. 

• Decays such as — > + neutral particles (e.g. — >■ nvr^, — * Avr", 
K±/vr±^/i±i/^). 

To reconstruct these tracks, the microvertex tracking algorithm initially looks for triplets 
of i?0 hits in the VD (there must be at least one hit in each layer), not associated to 
tracks reconstructed by the standard tracking algorithm. This allows a determination of 
the track parameters in the i?0 plane. A second step searches for unassociated z hits 
in the modules containing the R(j) hits. If there are at least two z hits for a R(j) triplet, 
the full set of parameters is given for that track. To improve the momentum resolution 
and to remove products of hadronic interactions in the beampipe, the results from the 
VD tracking were refitted forcing the track to originate from the primary vertex. The 
probability of this fit had to be greater than 0.1%. 

2.1.2 Efficiency correction procedure 

To be as independent as possible from the detector simulation, especially from the 
modelling of the VD, the efficiency for reconstructing a VD track has been deduced 
directly from the data as follows. 

Hadronic interactions in the detector material with at least two outgoing tracks were 
reconstructed by fitting the candidate tracks to a common vertex. The algorithm allows 
for an arbitrary number and charge configuration of tracks associated to the vertex^ The 
distribution of the measured positions of these vertices shows the material distribution of 
the Delphi detector (Figure |1]). If the incoming particle which caused the interaction is 
charged^, there is the possibility to reconstruct it as a VD track. To find these particles, 
the VD tracks were extrapolated to the radius of the interaction vertex and linked as 
the incoming track to the vertex if the difference in the azimuth angle between the 
reconstructed vertex and the VD track was below 1.7°. This cut was chosen to achieve an 
efficiency and purity for the linking close to unity. These linked VD tracks have not been 
used as candidates for E hyperons to suppress background from hadronic interactions. 
Finally, hadronic vertices within the polar angle acceptance of the VD were selected and 
the number of interactions with a link to a VD track was compared to the number of all 

^ The charge of the outgoing tracks may not sum up to the charge of the incoming track for hadronic interactions because 
an atomic nucleus is involved in the reaction. 

*This is most often the case since the production rate of neutral particles causing hadronic interactions (mainly K^, n 
and A") is small compared to the production rate of charged pions, kaons and protons. 
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Figure 1: Reconstructed hadronic interactions in the material of the Delphi detector, 
used for the determination of the correction factor for the VD track efficiency. Left: Rz 
view; Right: Radial projection. 



selected interactions. The efficiency for reconstructing a VD track is then given through 

Phv/vd ' ^hv/vd /^<. 
Phv ' hv ' ^link ' Jc 

where Nhv is the number of hadronic vertices, Nhv/vd the number of hadronic vertices with 
a linked VD track, phv and Phv/vd the purities of these samples^, eunk the efficiency to 
link the incoming VD track to the vertex and fc the fraction of charged particles causing 
hadronic interactions in the detector material. Assuming that the value of fc is the same 
in real data and simulation and taking into account the 6 dependence by multiplying 
with the 6 distribution of E hyperons in simulated events, the following correction factor 
for the reconstruction efficiency was deduced: 

r = ^ = 0.91 ±0.06, (3) 

VD 

where the error is systematic and comes mainly from the uncertainty in the fraction of 
charged particles causing hadronic interactions, which has been estimated in the following 
manner. In the region considered for the positions of reconstructed vertices the fraction 
of charged particles causing the interactions in the simulation was fc = 0.84. As a 
conservative choice, the error on this fraction was chosen to cover the range up to unity 
within three standard deviations, thus fc = 0.84±0.053. It is important to stress that this 
method does not rely on a precise modelling of the material distribution of the detector 
because the hadronic vertices were only used as candidate endpoints for the VD tracks. 



2.1.3 Reconstruction of the decay — > n7r* 

The candidates for the outgoing pion had to fulfil the following criteria: 

• Measured track length > 30 cm. 

• Ap/p < 1. 

^Phv a-nd Phv/vd were estimated directly from the radial distribution of the vertices shown in Figure 



5 



• IP jo IP > 3 in R(\) and where IP denotes the impact parameter with respect to 
the primary vertex and cx/p its error. 

• No associated VD hits in the two innermost layers. 

• The track must not originate from a reconstructed hadronic interaction. 

The z information of the VD tracks with at least two VD hits was not used in this 
analysis, because of differences in the association of two z hits between real data and 
simulation. Thus all VD tracks were treated in exactly the same way. 

To find the decay vertex of the S^, the intersections (normally two) in the i?0 plane 
of the VD track with the previously selected pion candidates were reconstructed. To 
select the correct intersection point, consistency of this point with the incoming and 
outgoing track was required which included removing intersections in the hemisphere 
opposite to the tracks. If both points fulfilled all these cuts, the one with the lower 
radius was chosen. The z coordinate of this candidate decay vertex was given by the 
z coordinate of the outgoing pion candidate at this -R0 position. The polar angle Q 
of the incoming track was calculated using the z coordinates of this vertex and of the 
primary vertex completing the track parameters of the VD track. The momentum of the 
(unreconstructed) neutron candidate could be computed allowing the calculation of the 
invariant mass of the S^-candidate. The combination was rejected if the calculated Q 
did not lie within the polar angle range covered by the modules of the VD defining the 
VD track. Both tracks defining the vertex must have the same charge. In order to reject 
background, the following additional cuts were applied: 

• The probability of the particle decaying within the measured flight distance, calcu- 
lated under the S~ hypothesis, had to be lower than 97%. This efficiently removes 
background from the decays K/vr — > /iz/^. Due to their long lifetime (cr(K^) ^ 
3.7 m and cr(7r^) ^ 7.8 m) charged kaons and pions decaying inside the sensitive 
volume of the tracking chambers tend to have very low momenta. Thus, their decay 
probability as defined above is close to 100%. 

• I cos 6'* I < 0.8, where d* denotes the angle between the outgoing track and the VD 
track, calculated in the rest frame of the particle reconstructed as VD track, where 
the pion mass was assumed for the measured outgoing track and the neutron mass 
for the undetected particle. 

The resulting mass spectrum of the E candidates is shown in Figure ^. A fit with 
Gaussians for the signal and the refiection from S~ — > Avr" (cr = 4.91 cm), which has 
the same signature of a kink|^, and a function of the form 

F[M) = ai((M - 1.079)(a4 - M))"^ exp(-a3(M - 0.5(1.079 + a^))) (4) 

for the background, where 1.079 is the sum of the masses of the neutron and the pion (in 
GeV/c^) and thus the lower kinematical limit, gives a mean mass of 1196.5 ±0.4 MeV/c^ 
and a width a = 12.3 ± 0.4 MeV/c^, in good agreement with the expectation from 
simulation of 1197.4 ± 0.3 MeV/c^ and 12.2 ± 0.6 MeV/c^ respectively. To subtract the 
background, two Gaussians with all parameters left free have been used for the signal to 
take into account the momentum dependence of its width. This results in a measured 
signal of 4820 ± 109 S decays and 870 ± 95 H~ decays (statistical errors only). 



®The decays f2 AK ,r}'n also show this signature, but their contribution is negligible due the low f2 production 
rate. 
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Figure 2: The invariant mass spectrum for the S candidates selected as described in 
the text. Dots are the data. The curve shows the result of a fit with two Gaussians for 
the signal, one Gaussian for the reflection from the decay S~ —>■ Avr" (dark shading) and 
a function of the form (^) for the smoothly varying background (light shading). The 
per degree of freedom of the fit is || = 1.24. 



2.1.4 Measurement of the S differential production rate 

Due to the larger branching ratio into the final state nvr (see Eq. (|l])) and the longer 
lifetime (crs- ~ 1.8 • crs+), the efficiency to reconstruct the decay S ^ nvr with the 
method described above is much higher for the S~ than for the To obtain the S~ 
cross-section 

e a 

where N denotes the number of signal events, e~ the efficiency for S~ reconstruction, 
a = ^ the ratio of the efficiencies for detecting a S~ or E+, and N-^+ the true number 
of S"*" in the signal, an estimate for the S"*" rate to be subtracted has to be made. For 
this, the ratio of to S~ production from the simulation (Jetset 7.4 with default 
parameters) has been assumed. An error of 20% has been assigned to this ratio. In the 
Xp bins {xp = 2p/^/s) where the deviation of this ratio from unity is greater than 20%, 
this deviation has been taken as the systematic error. 

The differential S~ production rate has then been measured in nine Xp bins. For each 
Xp bin the efficiency has been estimated using the simulated sample with full detector 
simulation, taking into account the correction factor of equation (^ which was assumed 
to be independent of Xp. To obtain the number of signal events, the mass spectra have 
been fitted using a Gaussian function for the signal, a Gaussian function with position and 
width fixed to the values obtained from the simulation for the reflection S~ — >■ Ayr", and 
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Figure 3: The mass spectra for the different Xp bins. Dots are the data, the sohd 
hne shows the result from the fit and the shaded histogram the background used for 
subtraction. 



a function with four free parameters of the form (^) for the smoothly varying background 
(see Figure The reflection has not been fitted separately for the last two Xp bins 
due to the large width of the signal. In the last Xp bin the signal has not been fitted, 
but the background obtained from the simulation (including the contribution from the 
S-refiection and normalised to data statistics) has been subtracted to get the number of 
signal events. 



2.2 Consistency checks and systematic errors 



In addition to the systematic error already mentioned in Section |2.1.2 



some more 



systematic checks have been done to test the consistency between real data and simu- 
lation. The quantities used for the selection of the candidates have been compared and 
good agreement was found. The mass spectrum has been fitted separately for negative 
and positively charged S candidates. One obtains good agreement within the statistical 
errors (2351 ± 73 and 2298 ± 78 signal events, respectively). Since a = ^ (the values 



8 



are given in Table |^) is a function of the decay length of S~ and S"^, the flight distance 
distributions for real data and the simulation have been compared for each Xp bin and 
consistency within the statistical errors has been found. 

Since this analysis uses tracks which have been reconstructed using only the microver- 
tex detector of Delphi, a good internal and external alignment of the VD is essential. 
It has been checked that the widths obtained for the signal in the different Xp bins show 
good agreement between real data and simulation within errors. The signal has been 
fitted separately for both z hemispheres defined through cos 9y£)_track ^ 0. Taking into 
account the 6 dependence of the efficiency for reconstructing a VD track (according to 
(I)), one obtains iVcos6»<o/^cos6»>o = 1-05 ± 0.05 (stat.). Another quantity sensitive to 
the alignment and effects of the tracking, such as misassociating VD hits to the outgoing 
pion, is the kink angle between the two tracks. The simulation describes the data well, 
even for very low values (< 5°) of this angle, where the reconstruction of two distinct 
tracks and the kink between them is most critical. 

To estimate the uncertainty of the number of signal events, the parametrisation of the 
background and the mass window for the background subtraction have been varied. 

The systematic errors due to the efficiency correction and the subtraction of the S"*" 
rate have been added linearly for the different Xp bins and thus treated as fully correlated. 
Statistical and other systematic errors have been added quadratically from bin to bin. 
Different systematic errors have been added quadratically. The contributions from the 
different sources to the total error are listed in Table |1|. 



error source 


absolute unc. 


relative unc. 


data statistics 


0.0016 


2.0% 


simulation statistics 


0.0012 


1.5% 


efficiency correction (equation (^) 


0.0043 


5.3% 


fit procedure 


0.0020 


2.5% 


S+ rate from simulation 


0.0028 


3.5% 


extrapolation to unobserved Xp region 


0.0081 


10.0% 



Table 1: The uncertainties of the E production rate. 



2.3 A(1520) Reconstruction 

Beyond the general cuts given in Section ^ for this analysis it was required that the 
track impact parameter to the primary vertex was less than 0.5 mm in the i?(^-plane and 
1 mm in the 2;-direction. This requirement strongly reduces contributions of tracks from 
particle re-interactions inside the detector material. Furthermore there must be at least 
two tracks inside the angular acceptance 47° < 6 < 133° of the BRICH. 

The pK~(pK"'") mass-spectra were then constructed for each bin of each individual 
kinematical variable using identified particles. Particle identification was performed com- 
bining dE/dx and BRICH information. According to the quality of particle identification 
the tagging categories loose, standard and tight tags are distinguished for each particle 
species as well as for so called "heavy" particles combining protons and kaons. To fur- 
ther improve the quality of particle identification for a track of given momentum and 
(assumed) particle type it was required that information from the detectors specified in 
Table (^ was present. 



9 





momentum range in GeV/c 


0.3 - 0.7 


0.7-0.9 


0.9 - 1.3 


1.3 - 2.7 


2.7-9.0 9.0-16.0 


16.0 - 45.0 


vr 


TPC 


LRICH S 


GRICH S 


K 


TPC 


LRICH S 


GRICH V 

+ 

LRICH S 


GRICH S 


P 


TPC 


TPC 

+ 

LRICH V 


LRICH S 


GRICH V 

+ 

LRICH S 


GRICH V 


GRICH S 



Table 2: Momentum ranges for particle identification: TPC denotes identification using 
the dE/dx measurement of the TPC, LRICH S (V) denotes identification using a signal 
(veto) of the liquid RICH, and correspondingly GRICH for the gas RICH. 



A particle was then taken to be a proton if it was tightly tagged. Kaons were required 
to be tightly tagged in the momentum ranges p < 3.5 GeV/c and p > 9.5 GeV/c. In the 
intermediate momentum range kaons were also identified by a tight heavy particle tag 
combined with at least a standard kaon tag. To suppress combinatorial background it was 
required that the kaon momentum was between 28% and 100% of the proton momentum. 
This condition has been left out for the determination of the A(1520) spin alignment. 



Figure (^) shows the pK~ mass spectrum of the overall dataset and Figure (^) for 
the scaled momentum range Xp > 0.4. A clear A(1520) signal is observed in both mass 
spectra at about the expected mass. The signal to noise ratio improves for the higher 
Xp range indicating that a proper measurement can even be performed for high A(1520) 
momenta. It should, however, be noted that here the A(1520) signal became poorly 
visible if the particle identification requirements were relaxed. It has been checked that 
there are no prominent reflections from known particle decays in the pK~ mass spectrum. 

Particle identification inefficiencies, detector imperfections and the different kinemati- 
cal cuts imposed for charged particle and event selection, were accounted for by applying 
the approach first described in |l^, developed in |lT6| , [T7| , p!8| , p!9| and outlined in brief below. 



In the present analysis, the mass spectra were described by an anticipated distribution 
function, /(M, a), of the invariant mass M. The parameters a were determined by a least 
squares fit of the function to the data. The function /(M, a) had three components: 



/(M, a) = f^{M, a) + ^ (M, a) + /«(M, a), (5) 

corresponding to the signal (S), background (B), and refiection (R) contributions respec- 
tively. 

The signal function, f^{M,a), described the resonance signal in the corresponding 
invariant mass distributions. For the pK~ mass distributions it had the form 

f{M,a) = aiPS{M)-BW{M,a2,a3), (6) 

where the relativistic Breit-Wigner function BW for the A(1520) is multiplied by the 
function PS{M) to account for the distortion of the resonance Breit-Wigner shape by 
phase space effects (see []15| for details). 
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Figure 4: a) Differential pK mass spectra for the overall measured energy range, b) 
pK~ mass spectrum for Xp > 0.4. The histograms represent the fit described in the text. 



The background term, f^{M, a), was taken to be of the form 

/^(M, a) = BGj,tset{M) ■ P{M, a), (7) 

where BGjetset{M) represented the background shape generated by JETSET which de- 
scribes the gross features of the real background and P(M, a) = 1 + a^M + a^M^ + 
QqM^ + a-iM"^ was a polynomial of order four (sometimes aj has been fixed to zero) in- 
troduced to account for possible deviations of BGjetset{M) from the real background. 
All pairs of charged particles which do not come from the resonance considered and re- 
flections in the invariant mass spectra were included in the definition of BGjetset{M). 
This parametrisation of the background was different from the analytical form used in a 
previous DELPHI analysis [|l5|,|l6|,|l7l|ll • 

The third term, f^{M, a), represented the sum of all the reflection functions (RFi): 

f\M,a)=Y,a,+jRFi{M). (8) 

i 

Reflections arising from particle misidentification and contributing to Equation 
were considered, for example when resonances in the vr+K", K^K~ and pvr" systems 
(K*°, $, A°) distort the pK" mass spectra. Due to the efficient particle identification of 
the combined RICH and TPC tags and to the high identification purity provided by the 
tight cuts, the influence of reflections of this type was found to be much smaller than 
without particle identification. 

The functions RFi{M) in Equation (H) were determined from events generated accord- 
ing to the JETSET model. The contributions of the reflections to the raw mass spectra 
defined by the function iV'^(a) (see Equation (|) ) were then obtained by passing these 
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events through the detector simulation. This also took proper account of the influence of 
particle misidentification. 

In each mass bin, M, the number of entries Nm{o) predicted by the function /(M, a), 
representing a sum of contributions from the resonance signal, background and reflections 
(see [0), is given by 

Ng{n)=CgY.Sg^A^J^{a\ (9) 

n 



/n(«)= / r{M,a)dM, (10) 

where G = S , B or R, and M„ is the lower edge of the n-th histogram bin in the distribu- 
tion of the variable M. The coefficients An characterise the detector acceptance and the 
losses of particles due to the selection criteria imposed, and the Cm take into account the 
contamination of the sample by particles from decays, wrongly associated charged par- 
ticles, secondary interactions, etc. The smearing matrix Smn represents the experimental 
resolution. The An, Cm and Smn were estimated separately for the resonance signal, 
background and reflection contributions using the detector simulation program. Due to 
differences in the detector performance and data processing in different running periods, 
the simulated events generated for these periods were taken with weights corresponding 
to the relative number of events in the real data. 

The best values for a were then determined by a least squares fit of the predictions of 
Equation (H) to the measured values, Nm, by minimising the function 

X' = Y1 - ^m{d)r/al + ^(a, - a,)V(Aa,)^ (H) 

m i 

where cr^ = Nm + cr^(iVm) and <7{Nm) is the error on Nm due to the finite statistics 
of the simulation used to evaluate An, Cm and Smn- The second sum in Equation (|TT1) 
constrains some of the fitted parameters to the values di it Acij taken from external 
sources, such as the normalisation of the refiection functions to the particle production 
rates taken from this and other LEP experiments, and the mass and width taken from 



the PDG tables |21[. The errors obtained from the fits thus include the corresponding 
systematic components. As a cross-check the A(1520) mass was also left free in the fit. 
This lead to a mass of 1517.5 ± 1.7 MeV/c^ fully consistent with the PDG value. 

In order to determine the full experimental systematic error of the A(1520) cross- 
section the following sources of uncertainty were considered. The infiuence of an imper- 
fect description of stable hadron spectra by the detector simulation was estimated by 
varying the charged particle selections leading to an associated systematic error of 3%. 
An error of 9% was assigned due to the imperfect description of A(1520) production by 
the fragmentation model. It was estimated by varying kinematical distributions of the 
A(1520) like the decay angle and momentum distribution. An uncertainty of 3% on the 
resonance cross-sections is due to the imperfect description of the particle identification 
efficiency and the error due to the branching ratios assumed is 2.2% [^]. Uncertain- 
ties due to the unknown momentum dependence of the refiections (3%) were assessed by 
reweighting the shape of the momentum distribution predicted by the model in a range 
allowed by the fit. Finally the uncertainty in the A(1520) line-shape, the background 
parametrisations and the infiuence of the binning was estimated to be 12% by choosing 
different parametrisations and by changing the binnings of the mass spectra. Adding the 
individual contributions in quadrature leads to a final relative experimental error of the 
A(1520) cross-section of 16%. 
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3 Results 
3.1 E- 

The results for the differential production rate and reconstruction efficiencies are 
listed in Table |[ The differential S~ production rate is plotted in Figure || together 
with the predictions of the Jetset 7.4 (default parameters) generator and a previous 
measurement of the OPAL collaboration The higher S"*" rate compared to the S~ 
rate in Jetset is due to secondary E hyperons, mainly from decays of charm particles. 
The shape of the Xp spectrum is well described in the simulation. 



Xp 


e-(%) 


a = f+ 




N Signal 




1 d<^-s- 




n.d.f. 




0.03 - 


0.06 


2.8 ± 


0.1 


6.5 


1.7 


860 


± 52 


0.671 ± 


0.041 ± 


0.064 


0.06 - 


0.08 


5.6 ± 


0.2 


4.5 


1.2 


745 


± 44 


0.422 ± 


0.025 ± 


0.042 


0.08 - 


0.10 


5.9 ± 


0.2 


3.4 


1.1 


725 


± 40 


0.361 ± 


0.020 ± 


0.038 


0.10 - 


0.12 


6.5 ± 


0.3 


3.4 


1.3 


560 


± 34 


0.252 ± 


0.015 ± 


0.032 


0.12 - 


0.15 


6.2 ± 


0.2 


3.2 


1.1 


630 


± 35 


0.195 ± 


0.011 ± 


0.026 


0.15 - 


0.19 


5.0 ± 


0.2 


2.7 


1.4 


470 


± 32 


0.128 ± 


0.009 ± 


0.020 


0.19 - 


0.25 


4.1 ± 


0.2 


2.7 


1.2 


450 


± 29 


0.097 ± 


0.006 ± 


0.016 


0.25 - 


0.35 


3.7 ± 


0.2 


2.2 


1.2 


380 


± 26 


0.051 ± 


0.004 ± 


0.010 


0.35 - 


0.50 


4.0 ± 


0.3 


2.2 




175 


± 17 


0.015 ± 


0.001 ± 


0.005 



Table 3: Efficiencies and differential S production rate in bins of Xp. The errors given 
for e~ and Nsignai are coming from the simulation and data statistics, respectively. 

Integrating the differential production rate in the measured Xp range from 0.03 to 0.5 
gives 

{N^-lNl""^) = 0.065 ± 0.002 (stat.) ± 0.006 (syst.) (Jetset 7.4: 0.054). 

To get the total production rate, the Jetset 7.4 simulated data has been used to 
extrapolate to the unobserved Xp range. In Jetset 7.4, 18% of the S~ hyperons are 
produced with Xp < 0.03 and 2% with Xp > 0.5. The production cross-section in this Xp 
range has been scaled by the ratio of the measured and simulated cross-section in the 
measured range 0.03 < Xp < 0.5. This scaled cross-section has then been added and a 
systematic error of 50% has been assigned to it. 

This gives the mean number of S~ hyperons produced in multihadronic decays 

{N^^/N^'"'^) = 0.081 ± 0.002 (stat.) ± 0.006 (syst.) ± 0.008 (extr.), 

where the last error comes from the extrapolation. 

This result is compatible with the corresponding rate in Jetset 7.4 {{N-£- / N^^'^) = 
0.068). It is also in good agreement with the measurement from Opal [TT[|: {Nj^~/N^^'^) = 
0.083 ± 0.006 (stat.) ± 0.009 (syst.). 

3.2 A(1520) 

The A(1520) has = |~ with isospin and quark content (uds). It decays strongly 
but with a comparably small decay width of 15.6 ± 1.0 MeV/c^ as it is a Dq^ state which 
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Figure 5: The measured differential production rate in comparison with the Jetset 
model and the OPAL measurement [jTT|. The data points are plotted at the mean Xp 
position in the corresponding bin. The statistical and systematic errors have been added 
quadratically. 



predominantly decays into a = | baryon and one pseudoscalar meson. The A(f520) 
branching fraction used in the analysis is taken as half of the branching fraction to nucleon 
and kaon of 45%: 

E(A(1520) ^pK") ^22.5% . 

The total A(1520) rate is measured from a fit to the mass spectrum given in Figure (^) 
corresponding to the scaled energyQ range 0.035 < xT"'^ < 1: 

(A^A(i520)/A^z'^) = 0.0285 ± 0.0048 (fit) . 

The fit error includes the statistical error and also accounts for uncertainties in the 
A(1520) mass and width and the normalisation of the reflection functions (see Section 
|2.3|) . The fitted rate agrees well with the integrated rate from the A(1520) inclusive xe 
spectrum (see below). To estimate the total rate of A(1520) production this value has 
to be corrected for the small unmeasured xe range < 0.035. This correction was 

7 A(1520) _ 2e*('-'S^°T 
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estimated using a modified Jetset 7.4 model (see below) normalised to the data in the 
measured momentum range. An additional systematic error of 50% of this correction has 
been assumed. This yields for the overall A(1520) production rate 



0.0293 ± 0.0049 (fit) ± 0.0047 (syst.) ± 0.0003 (extr.) 



This result is slightly higher 



(iVA(1520)/A^z'") 

including the systematic errors considered in Section 
than, though fully consistent with, the Opal result p. 

In order to exclude that the observed A(1520)'s originate predominantly from heavy 
(b,c) particle decays the A(1520) production has been determined individually for samples 
strongly enriched in b-quarks and light quarks respectively ||20|. This enrichment relies on 

No sig- 
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increased particle impact parameters due to the high lifetime of B hadrons 
nificant change in A(1520) production has been found in either sample which leads to the 
conclusion that the dominant part of A(1520) production originates from fragmentation. 

The observed A(1520) rate is comparable to that of the E*^(1385) [p|,p!0|,^ which has 
the same strangeness and total spin. It can be concluded that the orbital excitation {L = 
1) of the A(1520) does not lead to a suppression of particle production. The comparably 
large observed A(1520) rate suggests that also other orbitally excited baryonic states are 
produced in fragmentation. In consequence, due to the vast amount of orbitally excited 
states, a large part of the observed stable baryons may descend from these excited states. 
This agrees with the expectations of the phenomenological models 0J^. Thus for baryons 
the situation is likely to be similar to the mesonic case. It is highly desirable to verify 
the production of other orbitally excited baryonic states, however, this is experimentally 
demanding due to the large width and complicated decay modes of these states. 



A(1520) ^ p K" 



0.8 - 




-1 -0.5 0.5 1 

CO SI^K 



Figure 6: Differential A(1520) distribution as function of cos'^k for Xp > 0.07. The full 
line represents the result of a fit of the angular distribution, the dashed line represents 
the expectation for unaligned A(1520)'s and the dotted (dash-dotted) line represents the 
expectation for pi i + P_i_i = 1 (0). 
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In order to determine a possible spin alignment of the A(1520)'s the distribution of the 
cosine of the kaon angle, cos-i^i^, in the A(1520) rest system with respect to the A(1520) 
direction is plotted in Figure for Xp > 0.07. This distribution has been fitted with the 
expected form of the angular distribution: 

Ti// Q N l+ScOS^t^K , S-ScOS^t^K 

W{cosi^k) = P + ^ (12) 

with p = pii + p_i_i, Pa denoting the spin-density matrix element. The fit yields 
p = 0.4 ± 0.2, thus no significant A(1520) spin alignment is observed. The specified 
error includes the statistical error and the error due to the uncertainty in the background 
shape. 

Furthermore, the approximate symmetry of this distribution with respect to cos = 
presents evidence for the validity of this analysis as possible errors in particle identification 
and corrections would, in general, lead to a distorted cos{}k distribution. 




Figure 7: Inclusive A(1520) distribution as function of xe- Horizontal error bars indicate 
the bin width, vertical error bars are the fit errors described in the text. 

The measured scaled energy x^^^^"^^^ distribution is given in Table and shown in 
Figure (|^) compared to the Opal result 0. At low xe both measurements agree within 
error. For xe > 0.3 this measurement yields a rate about three times higher. Note that in 
this energy range a clear A(1520) signal is observed (see Figure (|^)). The measurements 
are compared to predictions of modified Herwig 5.9 and Jetset 7.4 models. A(1520) 
production has been implemented in these models either by replacing the S*°(1385) by 
the A(1520) in the case of Jetset 7.4 or by adding only the A(1520) to the particle hst 
in case of Herwig 5.9. The predicted rates should not be expected to be well reproduced 
by the models and the model predictions have been renormalised to the observed A(1520) 
rate. The general shape of the fragmentation function is reproduced well by both models. 
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a;£;-range 


{xe) 


1/iVe.t. 


J A7'A(1520) / J™ 

aly ^ ' /aXE 


X /Ndf 


0.035 - 


0.07 


0.052 


0.203 


± 


0.068 ± 0.032 


28 / 27 


0.07 - 


0.15 


0.108 


0.083 


± 


0.026 ± 0.013 


32 / 27 


0.15 - 


0.20 


0.174 


0.062 


± 


0.026 ± 0.010 


18 / 27 


0.20 - 


0.30 


0.247 


0.058 


± 


0.025 ± 0.009 


25 / 27 


0.30 - 


0.50 


0.390 


0.025 


± 


0.011 ± 0.004 


48 / 27 


0.50 - 


0.70 


0.590 


0.0044 


± 


0.0018 ± 0.0007 


32 / 27 



Table 4: Differential A(1520) distribution as function of tlie scaled A(1520) energy 
_^A(i52o)^ Tlie first error is the fit error, the second error the systematic error (see Section 
|2.3|) . The x^/^df specifies the quahty of the fits. 



In order to demonstrate the importance of A(1520) production in Figure (§) the ratio 
of A(1520) to A production is shown as function of the scaled momentum Xp. For this 



comparison the measurement of is taken as it covers a similar range in scaled mo- 
mentum like this A(1520) measurement. It is seen that at small Xp, A(1520) production 
is about a factor 20 less than A production. At large Xp, however, this reduces to a factor 
~2.5. Such a behaviour would be expected from general fragmentation dynamics due to 
the higher mass of the A(1520). An increase of this ratio is also expected if many A's 
stem from resonance decays. Finally it is interesting to note that the ratio of A(1520) to 
proton production is identical, within error, at low energies and in hadronic Z decays 
fas calculated from this result and 



3.3 Discussion on baryon multiplicities 

The measured and A(1520) production rates may now be more generally compared 
to the Lep average values for all baryons [^. In Figure (^a) the sum of the production 



rates of all states of an isomultiplet from the well known baryon octet and decuplet and 
for the orbitally excited isoscalar A(1520) are shown as a function of the corresponding 
particle mass squared, M^. In the case that not all states of an isomultiplet are measured 
at Lep, equal production rates for these states are assumed. Therefore the vertical axis 
of the Figure (^) is denoted by (21 + 1) (n), where (n) is the mean number of a given 
particle per hadronic Z° decay. It is seen that the mass dependence of the production 
rates is almost identical for the following sets of baryons: 

1. A^, A with strangeness S=0; 

2. E, E*, A and A(1520) with S=l and 

3. S, E* with S=2. 

Finally the Q~ rate is well predicted if the same mass dependence with an additional 
suppression for the higher strangeness (S=3) is assumed, as taken from the difference 
of the first and second or second and third set, respectively. All data points from these 
three sets and the Q~ are well fitted {x^ /ndf = 5.9/6) by the ansatz (see dotted lines in 
Figure (|a)) 

7-^(2/+ 1) (n) = A exp(-6M2). (13) 

The values of the fitted parameters are: A = 20.1 ± 1.5, 7 = 0.482 ± 0.022 and b = 
2.61 ±0.08 (GeV/c^)~^. As decays of high mass particles feed down to lower mass states. 
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Figure 8: Ratio of the differential Xp distributions of A(1520)'s and A's. A data are 
taken from |^ . 



it is to be expected that especially the slope parameter b for primary produced baryons 
differs from the fitted value. If the production rates are weighted by 7""^ a universal mass 
dependence is observed for all baryons (see Figure (Pd)). 

A similarly simple behaviour was found for scalar, vector and tensor mesons [^]. Note 
that for mesons the production rates per spin and isospin projection were used. This is 
implicit in the figures shown in p6[, as the experimental rates for mesons are customarily 
specified for each isospin state individually. The mass dependence for mesons, contrary 
to the baryonic case is exponential in mass, M, thus for mesons 



7" 



2J+ 1 



{n) = A exp(-6M) 



(14) 



where k is the number of s and s quarks in the meson. For baryons and mesons an 
almost identical value of 7 ~ 0.5 was found. A stronger fall-off with the mass for baryons, 
compared to mesons, is to be expected as baryons are produced in pairs. 



4 Conclusions 

The differential cross-section of the S~ hyperon has been measured in multihadronic 
Z decays by reconstructing the kink between the S~ and the outgoing pion in the decay 
nvr". The measured production rate is 

{Nj^-fN^"^) = 0.081 ± 0.002 (stat.) ± 0.006 (syst.) ± 0.008 (extr.) . 

This result is about 20% above the prediction of the Jetset 7.4 model, but nevertheless 
compatible. The shape of the differential cross-section is well described in this model. 
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Figure 9: a) Sum of the baryon production rates of all states of an isomultiplet as a 
function of the squared baryon mass. Data were taken from b) Sum of the baryon 
production rates of all states of an isomultiplet weighted by 7""^ as a function of the 
squared baryon mass. 

The differential cross-section of the orbitally excited A(1520) baryon has been mea- 
sured relying strongly on the particle identification capabilities of the Delphi detector. 
The A(1520) is the only measured orbitally excited baryon produced in fragmentation at 
Lep. Its production rate is found to be 

(A^A(i520)/iV|^'^) = 0.0293 ± 0.0049(fit) ± 0.0047(syst.) ± 0.0003(extr.) . 

This rate is similar to that of the S*^(1385) and suggests that also other orbitally excited 
baryonic states are produced in fragmentation. The shape of the A(1520) fragmentation 
function is well described by the Jetset 7.4 and Herwig 5.9 fragmentation model if 
A (1520) production is introduced. 

It has also been shown that the production rates of all baryonic states measured so 
far at Lep can be parametrised by a phenomenological law. 
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